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The Standard Model prediction for the magnetic moment of the muon requires a 
determination of the electromagnetic form factor of the pion at high precision. It 
is shown that the recent progress in tttt scattering allows us to obtain an accurate 
representation of this form factor on the basis of the data on e+e~ — > tt+tt". The 
same method also applies to the form factor of the weak vector current, where the 
data on the decay t — > ■iT~n'^ Vt are relevant. Unfortunately, however, the known 
sources of isospin breaking do not explain the difference between the two results. 
The discrepancy implies that the Standard Model prediction for the magnetic 
moment of the muon is currently subject to a large uncertainty. 

Talk given at the Workshop Continuous Advances in QCD 2002/ Arkady f est 
in honor of the 60th birthday of Arkady Vainshtein, Minneapolis, May 2002. 



1. Motivation: magnetic moment of the muon 

The fabulous precision reached in the measurement of the muon magnetic 
moment [1] aUows a thorough test of the Standard Model. The prediction 
that follows from the Dirac equation, fi = e h/2m^, only holds to leading 
order in the expansion in powers of the fine structure constant a. It is 
customary to write the correction in the form 

li=- (1 + a) . 1 

Schwinger was able to calculate the term of first order in a, which stems 
from the triangle graph in fig. la and is universal [2], 

a=^ + Oia') . (2) 

The contributions of 0{a^) can also unambiguously be calculated, except 
for the one from hadronic vacuum polarization, indicated by the graph in 
fig. le. It is analogous to the contributions generated by leptonic vacuum 
polarization in figs, lb, Ic and Id, but involves quarks and gluons instead 
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of leptons. All of these graphs may be viewed as arising from vacuum 
polarization in the photon propagator: 



gauge terms . 



g2{l + n(g2)} 

The expansion of the self energy function n(t) in powers of a starts with 

n(i) = an(o)(i) + a2n(i)(i) + ... 

It is normalized by 11(0) = 0. The leading term can be pictured as 



n(o)(t) = O + 



As shown in rcf. [3], the modification of the Schwinger formula (2) that is 
generated by vacuum polarization can be represented in compact form:'^ 



a / 1 — X 
a = — ax ■ 



2 2 

X m,, 



(3) 



i + Wx) i-x 

Expanding this formula in powers of a, we obtain 

pi 2 r1 

a = - / dx{l-x)-— dx{l-x) nW(t^) + 0{a^) 

Jo Jo 

The first term indeed reproduces the Schwinger formula (2), which corre- 
sponds to graph la. The term linear in XI^^^ [t) accounts for graphs lb to le. 



*The formula only makes sense in the framework of the perturbative expansion [4] . The 
contribution generated by an electron looj), for instance, grows logajrithmically at large 
momenta and tends to — oo in the spacclikc region, 

Hence 1+ct He"'' (t) contains a zero in the vicinity of t = — exp(37r/a) ~ — lO'''''' GcV'^ 
(graphs Ic, Id and le push the zero towards slightly smaller values). At academically 
high energies, the photon propagator thus develops a " Landau pole" , reflecting the fact 
that the U(l) factor of the Standard Model does not give rise to an asymptotically free 
gauge theory. In the present context, however, this phenomenon is not relevant - we are 
concerned with the low energy structure of the Standard Model. 



3 



The contribution involving the square of II^^^ (t) describes the one-particle 
reducible graphs with two bubbles, etc. 

The vacuum polarization due to a lepton loop is given by {£ = e, /j., r) 

T(o)..^_± r^,_^ 

37r Am? s{s-t) 



nr(i) = iz I d^^. > (4) 



Re{s) = \l ^ 1 + 



The hadronic contribution cannot be calculated analytically, but it can be 
expressed in terms of the cross section of the reaction e+e~ hadrons. 
More precisely, the leading term in the expansion of this cross section in 
powers of a, 

is relevant. In terms of this quantity the expression reads 

At low energies, where the final state necessarily consists of two pions, the 
cross section is given by the square of the electromagnetic form factor of 
the pion, 

^''(^) = i(l-^j 1^(^)1'' s<9M2. (6) 

Numerically, the contribution from hadronic vacuum polarization to the 
magnetic moment of the muon amounts to ai„p ~ 700 x 10^^''. This is a 
small fraction of the total, a = 11659 203(8) x 10~^° [1], but large com- 
pared to the experimental uncertainty: a determination of Ohvp to about 
1% is required for the precision of the Standard Model prediction to match 
the experimental one. Since the contribution from hadronic vacuum polar- 
ization is dominated by the one from the two pion states, this means that 
the pion form factor is needed to an accuracy of about half a percent. 
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2. Comparison of leptonic and hadronic contributions 

Graphically, the formula (6) amounts to 



TT 



There are three differences between the pionic loop integral and those be- 
longing to the lepton loops: 

• the masses are different 

• the spins are different 

• the pion is composite - the Standard Model leptons are elementary 

The compositcncss manifests itself in the occurrence of the form factor F{s), 
which generates an enhancement: at the p peak, |-F(s)p reaches values of 
order 45. The remaining difference in the expressions for the quantities 
Re{s) and Rh{s) in eqs. (4) and (6) originates in the fact that the leptons 
carry spin i, while the spin of the pion vanishes. Near threshold, the an- 
gular momentum barrier suppresses the function Rh{s) by three powers of 
momentum, while Re{s) is proportional to the first power. The suppression 
largely compensates the enhancement by the form factor - by far the most 
important property is the mass: in units of 10^^°, the contributions due 
to the e,/x and r loops are 59040.6, 846.4 and 4.2, respectively, to be com- 
pared with the 700 units from hadronic vacuum polarization. The latter is 
comparable to the one from the muon - in accordance with the fact that 
the masses of pion and muon are similar. 



3. Pion form factor 

In the following, I disregard the electromagnetic interaction - the discussion 
concerns the properties of the form factor in QCD. I draw from ongoing 
work carried out in collaboration with Irinel Caprini, Gilberto Colangelo, 
Simon Eidelman, Jiirg Gasser and Fred Jegerlehner. 

The systematic low energy analysis of the form factor based on chiral 
perturbation theory [5] has been worked out to two loops [6]. This ap- 
proach, however, only covers the threshold region. The range of validity 
of the representation can be extended to higher energies by means of dis- 
persive methods [7], which exploit the constraints imposed by analyticity 
and unitarity. Our approach is very similar to the one of de Troconiz and 
Yndurain [8]. For a thorough discussion of the mathematical framework, I 
refer to Heyn and Lang [9]. 
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We represent the form factor as a product of three functions that account 
for the prominent singularities in the low energy region: 

F(s) =G2(s) X G3(s) X G4(s) . (7) 

The index has to do with the number of pions that generate the relevant 
discontinuity: two in the case of G2, three for G3 and four or more for G4. 

The first term represents the familiar Omncs factor that describes the 
branch cut due to tt^tt" intermediate states (states with two neutral pions 
do not contribute, because the matrix element (tt'^tt'^ out | |0) vanishes, 
on account of Bosc statistics). The corresponding branch point singularity 
is of the type ImG2(s) ~ (s — S2) 2 , with = 4M^. The Watson final state 
interaction theorem implies that, in the elastic region, 4M^ < s < 9M^, 
the phase of the form factor is given by the P-wave phase shift of the elastic 
scattering process tt+tt" tt+tt". Denoting this phase shift by S{s), the 
explicit expression for the Omncs factor reads: 

^ / X \ s dx 6(x) 1 

\ J AMI X{X- S) ^ 

The function 63(5) contains the singularities generated by Stt interme- 
diate states: G3(s) is analytic except for a cut starting at S3 = 9M^, with 
a branch point singularity of the type Im G3(s) ~ (s — ss)^- If isospin sym- 
metry were exact, the form factor would not contain such singularities: in 
the limit m„ = rrid, the term G3(s) is equal to 1. Indeed, isospin is nearly 
conserved, but the occurrence of a narrow resonance with the proper quan- 
tum numbers strongly enhances the effects generated by isospin breaking: 
the form factor contains a pole close to the real axis, 

G3(s) = l + e^— + ... s^ = {M^-\iT^f . (9) 

S^jj s 

This implies that, in the vicinity of s = M^, the form factor rapidly varies, 
both in magnitude and in phase. The pole term cannot stand by itself be- 
cause it fails to be real in the spacelike region. We replace it by a dispersion 
integral with the proper behaviour at threshold, but this is inessential: in 

the experimental range, the representation for 63(5) that we are using can 
barely be distinguished from the pole approximation (9). 

Isospin breaking also affects the scattering amplitude. In particular, it 
gives rise to the inelastic reaction 27r Stt, with an amplitude proportional 
to m„ — md- Hence unitarity implies that, in the region 9AfJ < ,s < 16M^, 
the elasticities of the partial waves are less than 1. Numerically, the effect 
is tiny, however, because it is of second order in m„ — md- To a very 
high degree of accuracy, the first two terms in eq. (7) thus account for 
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all singularities below S4 = 16M^ - the function G4(s) is analytic in the 
plane cut from S4 to 00. Phase space strongly suppresses the strength 
of the corresponding branch point singularity: Im 6*4(5) cx (s — 54)2. A 
significant discontinuity due to inelastic channels only manifests itself for 

S>Sin = {M^ + M^)'^. 

We analyze the background term 04(3) by means of a conformal map- 
ping. The transformation 

Sin Si y/Sin S . , 

Z = . = , liUl 

VSin - Si + ^Sin - S 

maps the s-plane cut along s > Sin onto the unit disk in the z-plane, so that 
the Taylor series expansion in powers of z converges on the entire physical 
sheet, irrespective of the value of the arbitrary parameter si. We truncate 
this series after the first few terms, thus approximating the function G4{s) 
by a low order polynomial in z. 



4. Roy equations 

The crucial element in the above representation is the phase S(s). The main 
difference between our analysis and the one in ref. [8] concerns the input 
used to describe the behaviour of this phase. In fact, during the last two 
or three years, our understanding of the tttt scattering amplitude has made 
a quantum jump. As a result of theoretical work [10-12], the low energy 
behaviour of the S- and P-waves is now known to an amazing accuracy - to 
my knowledge, tttt scattering is the only field in strong interaction physics 
where theory is ahead of experiment. 

The method used to implement the requirements of analyticity, unitar- 
ity and crossing symmetry is by no means new. As shown by Roy more 
than 30 years ago [13], these properties of the scattering amplitude subject 
the partial waves to a set of coupled integral equations. These equations 
involve two subtraction constants, which may be identified with the two 
S-wave scattering lengths ag, a^. If these two constants are given, the Roy 
equations allow us to calculate the scattering amplitude in terms of the 
imaginary parts above the "matching point" Em = 0.8 GeV. The available 
experimental information suffices to evaluate the relevant dispersion inte- 
grals, to within small uncertainties [10,12]. In this sense, ag, Oq represent 
the essential parameters in low energy tttt scattering. 

As will be discussed in some detail in the next section, chiral sym- 
metry predicts the values of the two subtraction constants and thereby 
turns the Roy equations into a framework that fully determines the low 
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energy behaviour of the tttt scattering amphtude. In particular, the P- 
wave scattering length and effective range are predicted very accurately: 
a\ = 0.0379(5) Af-^ and h\ = 0.00567(13) M"! The manner in which 
the P-wave phase shift passes through 90° when the energy reaches the 
mass of the p is specified within the same framework, as well as the be- 
haviour of the two S-waves. The analysis reveals, for instance, that the 
isoscalar S-wavc contains a pole on the second sheet and the position can 
be calculated rather accurately: the pole occurs &t E = — \iT„, with 
M„ = 470 ± 30 MeV, = 590 ± 40 MeV [11], etc. 

Many papers based on alternative approaches can be found in the liter- 
ature. Padc approximants, for instance, continue to enjoy popularity and 
the ancient idea that the a pole represents the main feature in the isoscalar 
S-wave also found new adherents recently. Crude models such as these may 
be of interest in connection with other processes where the physics yet re- 
mains to be understood, but for the analysis of the nn scattering amplitude, 
they cannot compete with the systematic approach based on analyticity and 
chiral symmetry. In view of the precision required in the determination of 
the pion form factor, ad hoc models are of little use, because the theoretical 
uncertainties associated with these are too large. 



5. Prediction for the tttt scattering lengths 

Goldstone bosons of zero momentum do not interact: if the quark masses 

rriu, rud are turned off, the S-wavc scattering lengths disappear, aj], Oq — > 0. 

Like the mass of the pion, these quantities represent effects that arise from 

the breaking of the chiral symmetry generated by the quark masses. In 

fact, as shown by Weinberg [14], and Oq are proportional to the square 

of the pion mass 

7/W"2 M"^ 
" ™- +0{M^) , a^ = -^+0(M,^) . 



The corrections of order contain chiral logarithms. In the case of a^, 
the logarithm has an unusually large coefficient 



327rF2 1^ 2 (47rF^)2 

This is related to the fact that in the channel with 1 = 0, chiral symme- 
try predicts a strong, attractive, final state interaction. The scale Aq is 
determined by the coupling constants of the effective Lagrangian of O(p^): 

9 , Ag 20 40 X 5 5 
2^"Ai = 21^^+21^^-14^^ + 2^^+2 • 
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The same coupling constants also determine the first order correction in 
the low energy theorem for Og. 

The couplings £i and £2 control the momentum dependence of the scat- 
tering amplitude at first nonleading order. Using the Roy equations, these 
constants can be determined very accurately [11]. The terms £3 and £4, 
on the other hand, describe the dependence of the scattering amplitude on 
the quark masses since these cannot be varied experimentally, £3 and £4 
cannot be determined on the basis of nn phenomenology. The constant £3 
specifies the correction in the Gell-Mann-Oakes-Renner relation [15], 

M^ = M^\l--——^£3 + 0{M^)\ . (11) 



2 {AnFf 

Here stands for the term linear in the quark masses, 

= {mu + md)muu\0)\^ (12) 

{F and (Oj uu |0) are the values of the pion decay constant and the quark 
condensate in the chiral limit, respectively). The coupling constant £4 oc- 
curs in the analogous expansion for Fj^, 

A low energy theorem relates it to the scalar radius of the pion [16], 

The dispersive analysis of the scalar pion form factor in ref. [11] leads to 

{r\ = 0.61 ± 0.04 fm^ . (15) 
The constants £i,...£4 depend logarithmically on the quark masses: 

A? 

/. = In 't- i = 1 A 

H — in ^2 ! « — 

In this notation, the above value of the scalar radius amounts to 

A4 = 1.26 ±0.14 GeV . (16) 

Unfortunately, the constant £3 A3 is not known with comparable preci- 
sion. The crude estimate for £3 given in ref. [16] corresponds to 

0.2 GeV < A3 < 2 GcV . (17) 

It turns out, however, that the contributions from £3 are very small, so 
that the uncertainty in A3 does not strongly affect the predictions for the 
scattering lengths. This is shown in fig. 2, where the values of Oq, Oq 
predicted by ChPT are indicated as a small ellipse. 
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Figure 2. Prediction for the S-wave tttt scattering lengths 

6. Experimental test 

Stern and collaborators [17] pointed out that "Standard" ChPT relies on 
a hypothesis that calls for experimental test. Such a test has now been 
performed and I wish to briefly describe this development. 

The hypothesis in question is the assumption that the quark condensate 
represents the leading order parameter of the spontaneously broken chiral 
symmetry. More specifically, the standard analysis assumes that the term 
linear in the quark masses dominates the expansion of M^. According to 
the Gell-Mann-Oakes-Renner relation (12), this term is proportional to the 
quark condensate, which in QCD represents the order parameter of lowest 
dimension. The dynamics of the ground state is not well understood. The 
question raised by Stern et al. is whether, for one reason or the other, the 
quark condensate might turn out to be small, so that the Gell-Mann-Oakes- 
Renner formula would fail - the "correction" might be comparable to or 
even larger than the algebraically leading term. 

According to eq. (11), the correction is determined by the effective cou- 
pling constant £3. The estimate (17) implies that the correction amounts 
to at most 4% of the leading term, but this does not answer the question, 
because that estimate is based on the standard framework, where (0| uu |0) 
is assumed to represent the leading order parameter. If that estimate is 
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Figure 3. Interpretation of the data on the phase difference 5\ — in Generalized 
Chiral Perturbation Theory. 

discarded and is treated as a free parameter ("Generalized" ChPT), the 
scattering lengths cannot be predicted individually, but the low energy the- 
orem (14) implies that - up to corrections of next-to-next-to leading order 
- the combination 20,^ — Soq is determined by the scalar radius: 



The resulting correlation between and Oq is shown as a narrow strip 
in fig. 2 (the strip is slightly curved because the figure accounts for the 
corrections of next-to-next-to leading order). 

In view of the correlation between and Oq, the data taken by the 
E865-collaboration at Brookhaven [18] allow a significant test of the Gell- 
Mann-Oakes-Renner relation. The final state interaction theorem implies 
that the phase of the form factors relevant for the decay Tr+vr^e+Pe 
is determined by the elastic tttt scattering amplitude. Conversely, the phase 
difference 5q — 5\ can be measured in this decay. The analysis of the 4 • 10^ 
events of this type collected by E865 leads to the round data points in fig. 3, 
taken from ref. [19] (the triangles represent the Ke^ data collected in the 
seventies). The three bands show the result obtained in Generalized ChPT 
for ttg = 0.18,0.22,0.26, respectively. The width of the bands corresponds 
to the uncertainty in the prediction. A fit of the data that exploits the 
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correlation between afj and Oq yields 

= 0.216 ± 0.013 (stat) ± 0.004 (syst) ± 0.005 (th) [18] , 

where the third error bar accounts for the theoretical uncertainties. The 
result thus beautifully confirms the prediction of ChPT, = 0.220 ±0.005. 
The agreement implies that more than 94% of the pion mass originate in 
the quark condensate, thus confirming that the Gell-Mann-Oakes-Renner 
relation is approximately valid [19]. May Generalized ChPT rest in peace. 

7. Comparison of electromagnetic and weak form factors 

In the theoretical limit m« = and in the absence of the electromagnetic 
interaction, the vector current relevant for strangeness conserving semilep- 
tonic transitions is conserved. The matrix clement of this current that 
shows up in the decay r — » Tr~n° i^t is then determined by the electro- 
magnetic form factor of the pion. In reality, however, m„ differs from rud 
and the radiative corrections in r decay are different from those relevant 
for e~^e~ tt'^tt~. For the anomalous moment of the muon, r decays 
are of interest only to the extent that these isospin breaking effects are 
understood, so that the e.m. form factor can be calculated from the weak 
transition matrix element. 

The leading isospin breaking effects are indeed well understood: those 
enhanced by the small energy denominator associated with oj exchange, 
which are described by the factor Gs{s) introduced in section 3. As these do 
not show up in the weak transition matrix element, they must be corrected 
for when calculating the electromagnetic form factor from r decays. 

There is another effect that shows up in the process e+e~ — > tt+tt", but 
does not affect r decay: vacuum polarization in the photon propagator, as 
illustrated by the graphs below: 




The same graphs also show up in the magnetic moment of the muon: 
To avoid double counting, the data on the reaction e+e~ tt+tt" must 
be corrected for vacuum polarization, multiplying the cross section by the 
factor |1 + n(s)|^. In the timelike region, this factor is less than 1, so that 
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the correction reduces the magnitude of the form factor. In the oj region, the 

vaciiuni polarization due to Stt intermediate states generates a pronounced 
structure in n(s). While below that energy, the correction is of order 1%, 
it reaches about 7% immediately above the u) and then decreases to about 
3% towards the upper end of the range covered by the CMD2 data [20] . 

Unfortunately, applying the two corrections just discussed, the results 
for the form factor obtained with e^e~ collisions are systematically lower 
than those found in r decays [21-23]. The two phenomena mentioned 
above are not the only isospin breaking effects. Radiative corrections must 
be applied, to e+e" collisions [24] as well as to r decays [25] and terms of 
order m„ — need to be estimated as well. I do not know of a mechanism, 
however, that could give rise to an additional isospin breaking effect of the 
required order of magnitude. 

One way to quantify the discrepancy is to assume that the uncertainties 
in the overall normalization of some of the data are underestimated. Indeed, 
if the normalization of either the rate of the decay t tt^tt*^ or the 
cross section of the reaction e+e^ tt+tt" are treated as free parameters, 
the problem disappears. The renormalization, however, either lowers the 
ALEPH and CLEO data by about 4% or hfts the CMD2 data by this 
amount. 

While completing this manuscript, a comparison of the r and e"'"e~ data 
appeared [26], where the problem is discussed in detail. One way to put 
the discrepancy in evidence is to compare the observed rate of the decay 
r — > 'K~'K^Vr with the prediction that follows from the data on the reaction 
e+e~ — > 7r+7r~ if the known sources of isospin breaking are accounted for. 
Using the observed lifetime of the r, the prediction for the branching ratio 
of the channel r -> 7r-7r° Vr reads B^,,^ = 0.2408 ± 0.0031. The observed 
value, = 0.2546 ±0.0012, differs from this number at the 4(7 level [26]. 

A difference between the r and e+e~ results existed before, but with 
the new CMD2 data, where the hadronic part of vacuum polarization is 
now corrected for, the disagreement has become very serious. The problem 
also manifests itself in the values for the p width quoted by the Particle 
Data Group [27] : the value obtained by CMD2 is substantially lower than 
the results found by the ALEPH and CLEO collaborations. 
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So, unless the Standard Model fails here, either the experimental results 
for the electromagnetic form factor or those for the weak form factor must 
be incorrect. The preliminary data from KLOE appear to confirm the 
CMD2 results, but the uncertainties to be attached to that determination 
of the electromagnetic form factor yet remain to be analyzed. 

Currently, the discrepancy between the e+e" and r data prevents a test 
of the Standard Model prediction for the magnetic moment of the muon 
at an accuracy that would be comparable to the experimental value. The 
result for the contribution to the muon anomaly due to hadronic vacuum 
polarization depends on whether the e+e" data or the r data arc used 
to evaluate the electromagnetic form factor - according to rcf. [26], the 
corresponding central values differ by 17.2 x 10"^". Compared to this, the 
estimate (8 ± 4) x 10~^° for the contribution from hadronic light-by-light 
scattering [28-31] is a rather precise number. 

8. Asymptotic behaviour 

The behaviour of the form factor for large spacelike momenta can be pre- 
dicted on the basis of perturbative QCD [32]: 

The leading asymptotic term only involves the pion decay constant. The 

coefRcients B2, B4 . . . oi the fractional logarithmic corrections arc related to 
the pion distribution amplitude or null plane wave function ijj{x,fi), which 
is a function of the momentum fraction x of the quark and depends on the 
scale /U. Normalizing the wave function to the pion decay constant, the 
expansion in terms of Gegenbauer polynomials starts with 

{_ 50 3 _ 364 3 1 

Cl (z) = I (5^2 - 1) , C} {z) = f {21z^ - lAz-" + 1) , ... 

The wave function cannot be calculated within perturbative QCD and the 
phenomenological information about the size of the coefficients i32, . . . 
is meagre. It is therefore of interest to see whether the data on the form 
factor allow us to estimate these terms. 



The physics - in particular the sign - of the Ught-by-Ught contribution is well under- 
stood: The low energy expansion of this term is dominated by a logarithmic singularity 
with known residue [29,30]. In my opinion, the quoted error estimate is conservative. 
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In the representation (7), the asymptotic behaviour of the form factor 
can be accounted for as follows. One first continues the asymptotic formula 
(18) into the timclike region and reads off the asymptotic behaviour of the 
phase of the form factor: 

in A2 

If the asymptotic behaviour of the phase used for the Omnes factor agrees 
with this, then the Omnes formula (8) ensures that the ratio F{s)/G2{s) ap- 
proaches a constant for large spacelike momenta. The value of the constant 
is determined by and by the behaviour of the phase shift at nonasymp- 
totic energies. This implies that the background term G4{s) tends to a 
known constant for large values of s, or cquivalently, for z — 1. 

The corrections involving fractional logarithmic powers can also be ac- 
counted for with a suitable contribution to the phase. For the asymptotic 
expansion not to contain a term of order s~^/^, the derivative of G4(s) with 
respect to z must vanish at z — — 1. This then yields a representation of the 
form factor for which the asymptotic behaviour agrees with perturbative 
QCD, for any value of the coeffcients B2 and B4. 

Wc have analyzed the experimental information with a representation of 
this type, including data in the spacelike region, as well as those available 
at large timelike momenta [33-37]. The numerical analysis yet needs to 
be completed and compared with the results in the literature (for a recent 
review and references, sec for instance [38]). Our preliminary results are: 
If the fractional logarithmic powers are dropped {B2 = B4 = 0), we find 
that the asymptotic formula is reached only at academically high energies. 
With the value for B2 proposed by Chcrnyak and Zhitnitsky, the situation 
improves. For the asymptotic behaviour to set in early, an even larger value 
of B2 appears to be required. 

This indicates that the leading asymptotic term can dominate the be- 
haviour only for very high energies. A direct comparison of that term with 
the existing data, which only cover small values of s does therefore not 
appear to be meaningful. 

9. Zeros and sum rules 

Analyticity subjects the form factor to strong constraints. Concerning the 
asymptotic behaviour, I assume that | ln_F(s)| at most grows logarithmically 
for \s\ 00, in any direction of the complex s-plane. This amounts to the 
requirement that a) for a sufficiently large value of n, the quantity |F(s)/s"| 
remains bounded and b) the phase of F{s) at most grows logarithmically. 



15 



SO that the real and imaginary parts of -F(s) do not oscillate too rapidly at 
high energies. In view of asymptotic freedom, I take these properties for 
granted. If the form factor does not have zeros, the function 

F{s) 



In 



F(S2) 



(S2 - S)2 

is then analytic in the cut plane and tends to zero for \s\ oo. The branch 
point singularity at threshold is of the type ^{s) ~ (s — 52)"^ ■ Hence -0(5) 
obeys the unsubtracted dispersion relation 



*(«) = 



-/ 

71" Js- 



dx ■ , 



in the entire cut plane. The discontinuity across the cut is determined by 
the magnitude of the form factor: 



Im'^(s) = — - 



1 



In 



F{s) 



F{S2) 



S > S2 



Hence the above dispersion relation amounts to a representation of the form 
factor in terms of its magnitude in the timelike region: 



F(s) =|f (s2)|exp' 



(s2 - s)i 



f 



dx 



[x - S2) 2 (a; - s) 



In 



F{x) 



F{S2) 



(19) 



The relation implies, for instance, that the magnitude of the form factor in 
the timelike region also determines the charge radius.'^ 

Since the value at the origin is the charge, -F(O) = 1, the magnitude of 
the form factor must obey the sum rule 



ds 



Is2 S(S-S2)2 



In 



F{s) 



F{S2) 



= In \F{S2) 



(20) 



A second sum rule follows from the asymptotic properties. For the quantity 
In i^(s) not to grow more rapidly than the logarithm of s, the function 'tp{s) 
must tend to zero more rapidly than 1/s. Hence the magnitude of the form 
factor must obey the condition 



2M^ r 

71" Js 



ds 



S2 (5-52)2 



In 



F{s) 



F{S2) 



= 



(21) 



The relations (20) and (21) are necessary and sufficient for the existence of 
an analytic continuation of the boundary values of \F{s)\ on the cut that 



°For a detailed discussion of the interrelation between the behaviour in the spacelike and 
timelike regions, in particulaj: also in the presence of zeros, I refer to [39]. 
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(a) is free of zeros, (b) satisfies the condition F{0) = 1 and (c) behaves 
properly for large values of |s|. 

The above relations only hold if the form factor does not have zeros. 
In the scattering amplitude, zeros necessarily occur, as a consequence of 
chiral symmetry - indeed, the main low energy properties of the scattering 
amplitude may be viewed as consequences of the Adler zeros [40]. For 
the form factor, however, chiral perturbation theory implies that zeros can 
only occur outside the range where the low energy expansion holds: For 
the form factor to vanish, the higher order contributions must cancel the 
leading term of the chiral perturbation series. 

In quantum mechanics, the form factor represents the Fourier transform 
of the charge density. For the ground state of the hydrogen atom, for 
instance, the charge density of the electron cloud is proportional to the 
square of the wave function, which does decrease with distance, so that 
the corresponding form factor is positive in the spacelike region. It does 
not have any complex zeros, either. The wave functions of radially excited 
states, on the other hand, contain nodes, so that the form factor does exhibit 
zeros. Qualitatively, I expect the properties of the pion charge distribution 
to be similar to the one of the electron in the ground state of the hydrogen 
atom - in the null plane picture, the form factor again represents the Fourier 
transform of the square of the wave function [41]. In simple models such 
as those described in [42], the form factor is free of zeros. 

The hypothesis that the form factor does not contain zeros can be tested 
experimentally: The sum rules (20) and (21) can be evaluated with the 
data on the magnitude of the form factor. The evaluation confirms that 
the sum rules do hold within the experimental errors, but in the case of 
the slowly convergent sum rule (21), these are rather large. Alternatively, 
we may examine the properties of the form factor obtained by fitting the 
data with the representation (7). By construction, the first two factors in 
that representation are free from zeros, but the term 64(5) may or may 
not have zeros. In fact, as we are representing this term by a polynomial 
in the conformal variable z, it necessarily contains zeros in the z-plane 
their number is determined by the degree of the polynomial. The question 
is whether some of these occur on the physical sheet of the form factor, that 
is on the unit disk \z\ < 1. The answer is negative: we invariably find that 
all of the zeros are located outside the disk. It is clear that zeros at large 
values of )s) cannot be ruled out on the basis of experiment. In view of 
asymptotic freedom, however, I think that such zeros are excluded as well. 
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10. Conclusion 

The recent progress in our understanding of tttt scattering provides a solid 
basis for the low energy analysis of the pion form factor. The main prob- 
lem encountered in this framework is an experimental one: the data on 
the processes e^e~ — > tt+tt" and r tt^tt" Vt are not consistent with our 
understanding of isospin breaking. If the data are correct, then this repre- 
sents a very significant failure of the Standard Model - or at least of our 
understanding thereof. The discrepancy must be clarified also in order for 
the accuracy of the Standard Model prediction to become comparable with 
the fabulous precision at which the magnetic moment of the muon has been 
measured. 
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